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16. ASitiwi 


Memory effects in turbulent flow are studied. Based on 
experimental investigations of the wake flow of a hemi- 
sphere and cylinder, it is shown that such memory effects 
can be substantial and have a significant influence on 
momentum transport. Memory effects are described in terms 
of suitable memory functions. 
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MEMORY EFFECTS IN TURBULENCE 


J. 0. Hinze 

It Is shown that the effects in a turbulent flow are not 
strictly local and that memory effects play an essential role. 
Thus, momentum-transfer processes are not determined by the local 
gradient of the mean velocity alone. 

The Reynolds stresses can be thought of as consisting of two 
parts. The first part is the stress that would occur if it were 
determined solely by the local state of the flow, expressed in the 
local mean velocity gradient. The second part is an additional 
contribution from a memory of the influence of a change in this 
local gradient in the memory region of interest. 

A transfer of the gradient type will often occur in fully de- 
veloped flows, because either the additional memory effects are 
relatively small, or the flow is self-preserving. 

These additional memory effects can be directly detected in 
the development regions of turbulent flows, where the turbulence 
has no opportunity to adapt to new local flow conditions due to 
the rapid fluctuation of the main flow direction. 

The memory effect can be described with a suitable memory 
function . As an approximation, a simple exponential function is 
assumed for this. In this case the Reynolds stress is subject to 
a relaxation equation whose solution yields the axial distribution 
of the stress. 


* Numbers in the margin indicate pagination in the foreign text. 
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Based on experimental studies of the wake flow behind a hemi- 
sphere in a turbulent boundary layer as well as the wake flow of 
a cylinder, it is shown that such additional memory effects can 
be substantial. Both cases involve a flow in the development re- 
gion which fluctuates rapidly in the main flow direction. If the 
shear stress under consideration is expressed by the velocity 
gradient and an eddy viscosity, this eddy viscosity is equal to the 
value of the fully developed flow after correction for the addi- 
tional memory effect. 

1. General 


Flows which exhibit memory effects have been known for some 
time, mainly for non-Newtonian fluids. In this case the fluid 
should at least have viscoelastic properties which, as in rheo- 
logy, are determined by its molecular stvucture. 

But memory effects are also present in the turbulent flow of 
Newtonian fluids. The phenomenological theories of Prandtl and 
Taylor on the momentum transfer associated with turbulent or Rey- 
nolds shear stress essentially involve a memory effect determined 
by the Lagrangiar. autocorrelation of fluid particles in turbulent 
motion. Subsequent studies have shown that various phenomena in 
turbulence can be explained in principle by the assumption of a 
viscoelastic turbulence behavior , as in the formation of secondary 
flows in turbulent flow through a pipe of noncircular cross- 
section. It is apparent, then, that turbulence may be regarded as 
a hypothetical non-Newtonian fluid . A basic departure from rheo- 
logy must be made, however. 

First, the principle of local action is basic to rheology. 
According to this principle, the stresses in a particle are deter- 
mined or influenced by motions of the medium outside a certain very 
small material environment of the particle. "Particle" here re- 
fers to a . small quantity of matter with a specific identity and a 
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volume which Is very small In relation to the macroscopic length 
scale of the motions. The state of stress at a particular point 
In time Is determined entirely by the history of the motion In 
the Immediate environment of the particle up to that time. 

For " rheological ly simple * 1 fluids, the deformation history is 
described completely by the local deformation gradient (velocity 
gradient). Simple fluids are Isotropic, and they "remember" the 
past only via the deformation tensor up to the point of time under 
consideration. 

Second, a constitutive equation , i.e., a relation between 
deformation and stress, must satisfy various conditions. The 
equation must be Invariant with respect to a transformation of 
the coordinate system . It should also be Invariant with respect 
to time-dependent rigid motions such as translation and rotation ; 
In other words. It should be objective . 

A century ago, Bousslnesq regarded turbulent flow as a 
hypothetical fluid, specifically a Newtonian fluid. He assumed 
a linear relation between the stress tensor and the deformation . 
tensor by using an effective "eddy viscosity" as a scalar quantity 
With a small extension, Boussinesq's relation can be written 


— «/*< 


-i/VA„ + 
e 



(l) 


where P t 3 9 u j u i is the mechanical pressure, 
kinematic eddy viscosity. 


and e is the 
m 


Because the turbulent velocity fluctuations u^ are differ- 
ence velocities and thus objective, as Is the deformation tensor, 
relation (1) is objective. 
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If turbulence Is regarded as a Theologically simple fluid with 
non-Newtonian properties, relation (1) can be further supplemented 
as follows: 


**(*», t, V t , P,lg, dU„IQx m ) . ( 2 ) 

It can be shown that if the objectivity conditions are met, 
the following relation can be derived [1]: 


- UfU, — ej„ + 2 e m D„ + 4 e^D^Dn 


( 3 ) 


where Dfi „ 1/2 (e?7,/8*, + WJtx,). 

The quantities e n , e and e are still functions of the 

u _m . c 

principal invariants of and of the possible scalar invariants 
such as turbulent pressure. This nonlinear relation is analogous 
to the rheological equation of state for viscous fluids of Stokes 
and of Relner-Rlvlln . This relation satisfies the condition of 
local action but contains no "memory . " 

For a local action in turbulent flows, one must also consider 
very small particles which do not exceed the Kolmogroff micro 
length scale in order of magnitude. Now it happens that such 
small particles lose their identity (momentum) too quickly to 
make a significant contribution to turbulent transport. To be 
effective in this regard, the lumps of fluid must be rather large 
— of the same order of magnitude as the energy-containing 
eddies, i.e., the order of magnitude of the integral length scale. 
The action is essentially no longer local, and so higher deriv- 
atives of the mean velocity must also be considered. The radius 
of action is taken to be proportional to the eddy size under 
consideration. In this case the size of a unit area for momentum 
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transport (UjU^) can no longer be taken as small, but must be of 
the order of magnitude of these lumps of fluid. 

Memory effects In turbulence , expressed In a memory or relax- 
ation distance, are determined in part by the size of the eddies . 

A classic example of the dependence of the relaxation or memory 
distance on eddy size is found in the studies of Clauser [2]. 

Pig. 1 shows the result of these studies. A rod (d « 12.5 mm) 
is placed in a turbulent boundary layer of thickness 6 ■ 235 mm 
in a cross-stream direction, thereby producing a flow which is 
characterized by the local velocity defect AU^. The decay of this 
defect occurs much more rapidly near the wall (x 2 /6 “ 0.16; 
smaller eddies) than at a point more distant from the wall 
( x 2 /6 ■ 0.59; larger eddies). 




Fig. 1. Decay of a mean- 
velocity disturbance in a 
turbulent boundary layer. 


Fig. 2. Spectral energy distribu- 
tion of turbulence along a pipe axis. 

Key: a - Relative constriction 0.25; 

b- Without constriction; c - With 
constriction . 
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In his wave theory of turbulence, Landahl [3] showed that 
the only suitable length scale for the spatial behavior of an eddy 
Is the size of the eddy Itself . A wave component loses its identi- 
ty at a distance of approximately six of its own wavelengths. 
Measurements in a boundary layer showed that at X 2/6 * 0.5, the 
Lagrangian integral length scale is approximately 4 6 . Thus, the 
disturbance will decay completely only after about 24 6 . The de- 
pendence of the decay of velocity disturbances on eddy size was 
nicely confirmed by the experimental results of Lissenburg [4], 

The turbulent flow of air in a straight pipe of circular cross- 
section (D ■ 20 mm) was disturbed by a local constriction of the 
pipe. At various distances downstream of the constriction, mea- 
surements were made of mean velocity distribution and the inten- 
sity of turbulent velocity along the pipe axis and close to the 
wall. The Reynolds number relative to diameter and mean velocity 
was approximately 5000. 

The results of this experiment are shown in Fig. 2. After 
x ■ 20 D downstream of the constriction, the spectral energy dis- 
tribution is already equal to the undisturbed flow for frequencies 
higher than 100 s”^. For lower frequencies, the spectral distri- 
bution is not yet equal to the undisturbed flow even after 40 D. 

If, as an approximation, the propagation velocity of the distur- 
bances is assumed to equal the mean velocity, Landahl* s estimate 
for 20 D (« 0.4 m) leads to a frequency of 70 s _1 , which is of the 
same order of magnitude as 100 s”^. 

It can be concluded from Landahl *s work that the more the 
transfer of size is governed by the fine structure of the turbu - 
lence, the more local is the effect . This validates the assump- 
tion of a transfer of gradient type which is made in many new tur- 
bulence theories in order to determine the influence of velocity 
correlations of third and higher order. 
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Because the larger eddies are effective In momentum transfer, 
with an action that Is not limited to small regions, it may be 
asked whether the requirements of objectivity, and especially the 
condition of invariance under rotation, should still be satisfied, 
because centrifugal and Coriolis effects can no longer be 
neglected. Therefore, Lumley [1] has suggested that a constitutive 
equation for turbulence as a hypothetical fluid no longer bo 
subjected to this last condition. 

For the momentum transfer u^, we may write the following 
general expression for the time-averaged transport through a unit 
of area at right angles to the x^-dlrectlon : 

T 


The value U^(t Q ) is determined by the prior history of the 
particle for all t <. t Q . This prior history is determined in 
part by the distribution of mean velocity in the entire region 
under consideration and by the scale of the interaction between 
the particle and its environment during their respective motions. 
Averaged over many particles, this is associated with the 
Lagrangian autocorrelation. In principle, this aftereffect can 
be expressed in terms of a suitable memory or relaxation function. 
But the determination of such a suitable function generally 
presents a difficulty which has not yet been overcome. 

For a particle which has been in motion for a time t and has 
traveled a distance y^, a series expansion of U^Ct^jt) yields 


» 0 - (Pa - » (fcl + J - 5 (e£^.l 


( 5 ) 


If the first derivative of the mean velocity is sufficient, rela- 
tion (4), with 




yields 


i i _ 

_ ^ . Q j it .,<«*.«. -"n - ^ ^ ( 6 > 

^2k^L " u 2^o^ u k^ t o” t ^ is the La & ran S lan autocorrelation, and 
is an eddy viscosity which can have various values for the 
various directions x^. 

For t ®, we obtain the relation of Boussinesq, but with 
various possible values of the eddy viscosity The relation 

is not objective . Since the relation contains only the first ve- 
locity derivative, it appears as if the action is local in charac- 
ter. However, it must be noted that the appearance of the La- 
grangian autocorrelation points to an influence of prior history. 
This Influence has been nicely illustrated by Philip [5]. He 
imagines the velocity of a particle at each point in time as con- 
sisting of two parts. The first part corresponds to that which 
the particle still "remembers" of its velocity u(0) at time t ■ 0, 
while the second part consists of the momentum contributed by the 
exchange with its momentary environment: 

u(t) ■ u(0) R^(t) + u*(t). 

R T (t) is the Lagrangian autocorrelation coefficient and describes 
the influence of the past. 

For a stationary homogenous turbulence , it should be true that 
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u 2 (t) ■ u T (0) ■ const. Philip calls the first part Uj(t) ■ 
u(0) R^(t) the drift velocity . For very large times t, the 
memory of the original velocity becomes vanishingly small. 

It is assumed in the simple relation (6) that the velocity 
gradient in the region through which a particle passes during 
the memory period is sufficiently constant. In most turbulent 
flows of interest, however, this is not the case. Corrsin [6] 
recently made a critical study of the conditions which must be 
satisfied so that the assumption of a turbulent diffusion of the 
gradient type may be made . It Is shown that of the various con- 
ditions, the most Important Is that the length scale of the 
transport mechanism should be small relative to the distance 
over which the curvature of 3lJ,/3x 2 , l.e. 3 , U 1 /3xj, undergoes a 
substantial change . Corrsin considered a one-dimensional diffu- 
sion process in the x 2 -direction and a spatially symmetrical 
transport mechanism, so that the term with the second derivative 
is absent In a series expansion such as eqn. (5). Experience 
has shown that the error in assuming a diffusion of the gradient 
type is not large. This ic due to the fact that the diffusion 
is essentially one-dimensional in most cases. 

In the cases of three-dimensional transport considered here, 
the terms with the second derivative must also be considered in 
series expansion (6) as a second approximation, without assuming 
a diffusion of the gradient type. If this is not done, a non- 
negligible error can be made locally in the distribution of the 
mean velocity. As early as 19^2, Prandtl [7] pointed this out in 
a modification of his mixing-length theory for the environment of 
the maximum velocity of a free .Jet, in which the first velocity 
derivative vanishes. In this modif lea 4- • on he proposed a statis- 
tical mean for the first and second derivatives. 

Based on his measurements of turbulent diffusion of the 
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kinetic energy of the three turbulent velocity components In the 
wake flow of a cylinder, Townsend [8] concluded that there are 
places in the wake where diffusion takes place against the gradi- 
ent — specifically, near the energy-distribution maximum of the 
turbulence in a cross-section of the wake flow. 

It is also known that near a mean-velocity maximum in the 
case of asymmetrical flows, the turbulent shear stress and the 
mean-velocity gradient may have opposite signs, which car. be 
interpreted as a ’’negative energy production” of the turbulence 
[9-13]* The explanations advanced by Beguier [12], Hinze [14], 
Eskinazl [15] and other authors are all based on a non-negliglble 
influence of large eddies with a memory effect that requires at 
least that the second derivative of the mean velocity be taken 
into account. 

Apparently the turbulent shear stress, or more generally, 
the Reynolds stresses, cannot be regarded as a local property of 
the flow. This led Philips [16] to suggest that one should regard 
the gradient of the Reynolds stress, rather than the Reynolds 
stress itself, as a local property. 

In the cases previously mentioned, one must take into 
account the second velocity derivative, l.e., the change in the 
velocity gradient perpendicular to the main flow. But there are 
also cases in which the change of this gradient in the main 
flow occurs relatively quickly with a time or length scale which 
is no longer large in comparison with the scale of the large 
energy-containing eddies . In this respect the turbulence is no 
longer ’’adapted” to or ”in equilibrium” with the local conditions 
of the flow. 

In his study of boundary layers which undergo rapid spatial 
changes in the main-flow direction, Deissler [17] showed that the 
Reynolds stress remains approximately constant along one streamline 
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and may be equated with the initial value. In these necessarily 
short boundary layers, the memory distance extends over the entire 
length of the boundary layer considered. 

Other examples in which memory effects of the type mentioned 
have been detected are the wake flow behind a hemisphere in a 
turbulent boundary layer and the wake flow immediately behind a 
cylinder. These two cases will be discussed in greater detail 
below. 

The memory effects can be allowed for if the series expansion 
of Uj according to eqn. (5) is continued in relation (4), and 
first and higher derivatives of are also considered. 

Relation (4) then becomes 




d t' |^(/ e ) «*(<, - t’) - j u,(/,) u k (t a _ t') y, ((,;(')— + 


( 7 ) 


8 * 




+ 6 «,«.) Ufc - n nur-Ct ; 0^5- + - j Tx ; • 


In mos' ~ases it is sufficient to take into account only the 
terms up to the second derivative in a series expansion as a 
second approximation . In a manner analogous to relation (6), we 
find from eqn. (7) 


, V 0^1 if 


Qtrj 

d*' Oyi(f 0 ;O 1 


Sx t dx t 
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The last term In an approximation based on the assumption 
that the integral of a triple correlation can be written as the 
product of a suitable length Lj(t) and an eddy viscosity. A 
further approximation is that a scalar eddy viscosity is suffi- 
cient for e 2 ^ and £°k* because it is also assumed that t -+ 


- «*«i 




Wi l 

Sins*,] 


( 8 ) 


This equation can be generalized to a relation for -u^u^, 
one which is symmetrical in j and i. Because STJ^/ax.^ = 0 for an 
incompressible medium, we write the following expression so that 
the relation for UjU^ contains no contradiction: 


_____ p r i a 1 

— UfU | = — ~ + 2 1 e m Df t — tmLi 


(9) 


The second velocity derivative in eqn. (8) and eqn. (9) may 
signify a nonlocal action . This depends on what the length scale 
Lj is determined by. At right angles to the mean-flow direction. 
Inclusion of this term signifies a nonlocal action. In the main- 
flow direction, however, Lj can become relatively large due to 
the high convection velocity, and the action may be local in 
character. 

The influence of a spatial change in D ^ i on momentum 
transfer can also be taken into account with a memory function. 
Instead of eqn. (9) we then consider 


- u,u, = - + 2 fdf e „(0 D n (t - n 0(t') . 

Q o 


(10) 
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For large times t, the normalized memory function G(t) decays to 
zero. It is characteristic of the turbulent flow under consider- 
ation. Eqn. (10) points essentially to a local action. As 
mentioned earlier, however, relatively large lumps of fluid play 
a role in momentum transfer. If varies noticeably in the 
range of the order of magnitude of these lumps, the action is 
not local. However, a description based on eqn. (10) may be 
possible if the transfer takes place as if the action were local. 

Relation (10) is similar to that for a viscoelastic fluid. 
Crow [ 1 8 ] has suggested such a relation for describing the proper- 
ties. of the fine structure of a turbulent flow. 

2. Theory 

In this section we shall examine relations (9) and (10) more 
closely. Of particular interest for the cases investigated 
experimentally is the shear stress per unit mass (-u 2 u l ) , 
because these cases involve nearly one-dimensional flows with a 
pronounced main-flow direction. The mean flow is two-dimensional 
with the velocity components Uj in the main-flow direction x lS 
and U 2 in the direction x 2 at right angles to the main-flow 
direction, where Uj >> U 2 . The mean flow is independent of the 
x 3 -direction. 

As a second approximation, it is assumed that the memory 
effect of SUj/aXj, i.e., the influence of a relatively rapid 
change of d\J 1 /dx 2 , is sufficient in the main-flow direction. 
According to eqn. (9), 



with Aj = 1/2 . 
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Analogously, eqn. (10) gives 
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- «,«,«) = J di' t m (n ^-(t- n o(n . . 


(12) 


The last expression is more general than eqn. (11). The 
significance of these relations is that the shear stress may be 
envisioned as consisting of two parts, hence: 




The first part (-u u ) is the shear stress which would - 

2 1 g 

occur if it were determined entirely by the "local" state of the 
flow. Here "local" does not necessarily refer to a small region. 
It may be said that this part is "in equilibrium" with the local 
state. It may, for example, be equated with the shear stress 
according to the Boussinesq relation e m if the memory 

remains confined to a region in which 3U x / 3 x 2 is practically 
constant. 


The second part (-u 2 Uj ) is an additional contribution which 
a lump of fluid "remembers" during the transfer process from the 
influence of a variable 3 Uj/9x 2 in the memory region considered. 


Usually the relative intensity of the turbulence is not 
large, and it may be assumed that Taylor’s hypothesis of a frozen 
turbulence is valid: x t ^ Ujt. Eqn. (12) can then be written 

in a Eulerian description: 




= J dx[ e m (x') ( x't ) Q(x l - x\) . 


(13) 
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This expression contains the contributions (-u 2 u, ) and (-u 2 Uj ) g . 

If e m 3U 1 /3x 2 - const. In the region In which G(Xj - xj) > 0, we 
obtain with - x[) = 1 the Bousslnesq relation, and (-u 2 u x ) =0. 

As mentioned, eqn. (13) Implies essentially a local action. 

In the present work this relation Is successfully applied to se- 
veral turbulent flows. In these cases It Is apparently correct to 
assume that the momentum transfer is local In character but has a 
memory. 

An obvious step Is to associate the memory function with a 
Lagranglan autocorrelation. Although in principle the Lagrangian 
autocorrelation is not equal to a simple exponential function, 
such an assumption is sometimes made with satisfactory results in 
order to simplify calculations. We shall make the same assumption 
here for the memory function; thus, in normalized form, we obtain 




(14) 


is 'then an effective memory length . With this expression for 
G, eqn. (13) becomes 




J e m (x[) — ^ (zj) exp ^ — 


*i -*i V 

A > 


(15) 


Owing to the exponential function, expression (15) can now be 
regarded as the solution of an inhomogenous differential equation 
of first order, which is compatible with the relaxation process 
[19] considered above. This differential equation is 


A (' ~ «I«l) + (- «.«l) = ~ . 


( 16 ) 
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Up to now we have assumed the length to be Independent of 
Xj, However, experiments have shown that A x undergoes noticeable 
changes over distances equal to Its length. Equation (16) can be 
expanded In the sense that A t may still be a function of x t . The 
solution is then 


- = - «a«i(*u) ex P 


[ _ /^5 + / di: 


, *V 8 ^1 / - 1 


*1 

-/ 


dzi 

A(*i) 


(17) 


The integrals in the exponential expressions can be evaluated 
if we assume a power function for h l : 


A(*i) = • 


( 18 ) 


Solution (17) for n = 1 then becomes 
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• i-wi!* • (!*) t ~ (!*)“» !*P / (if) * ( 8,f )!^ - = t 1 *) 1 **" - 

y?0 ./ »/i- v 

■j 

Foi n / 1 we obtain 


(19) 


r x \-n _ x} 0 “*l 

= - UtU^Xu) exp I - ‘ a(1 _ n) -J + 

+ r [ ■ - *r^o] • / dli ^ e,p [siT^j] 


( 20 ) 


An equation for the term 3/9Xj (-u 2 u ) , which occurs in eqn. 
(16), is also derived from the dynamic equation for the Reynolds 
stress, specifically f20]: 
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The second term on the right side Is usually negligibly small. 
Disregarding viscosity effects, the equation reduces to 


«.«■>=* * 5 ^- 



• r the pressure term , Rotta [21] suggested the following 

combination with the shear stress: 


4. = 

Q ' 0 *, 8 x,/ 



(- «,«,) , 


where q ? = (q 2 ) 1 ^ ■ (u^u.^) 1 ^ ; L is a length scale of the size of 
an integral length scale; and A is a numerical constant. Hence, 


8 u\ 8tf, A <f . 


( 21 ) 


L is, as mentioned, of the size of the integral length scale. For 
this we select A^, as an axial Integral length scale and set L = 
const A lf ,. We also set A x = U l I lL = (UJu'i) A lL = const £7,/uJ A v with uj = 
(u 2 ) 1 / 2 ; I^ L is the Lagrangian integral time scale, A^ L is the 
Lagrangian integral length scale. If now A/L q’/Uj = const* 
•q'/uJ’l/Aj, we obtain eqn. ( 1 6 ) from eqn. (21) if const q'/u| = 1.. 
Further let u 2 A./U, = u 2 I,, = e I, T /I OT , and for the Lagrangian 
integral time scale, 1^ 'v* ^IL* which is the case only in an iso- 
tropic state of turbulence and is true only in rough approximation 
for a nonisotropic turbulence. In any case, the relaxation 
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equation (16) is not incompatible with the dynamic equation for 
shear stress. 

It is our intention to show memory effects in flows which 
fluctuate rapidly in the mean flow direction . It may be asked 
whether a memory effect can be present in flows which do not 
fluctuate rapidly in this way. According to relations (11) and 
(16), the memory effects are negligible if 

< eC/./ftr, ( 22a ) 


or if 

4|8(- i < “ “i*! • ( 22b ) 

Now there is a class of flows for which these conditions are 
not satisfied, and yet there is no demonstrable memory effect. 
These are flows which possess the property of self-preservation . 

Self-Preserving Flows 

These are flows with a pronounced main flow direction, whose 
structure maintains its similarity during development or decay in 
this direction . Many free turbulent flows such as free jets and 
wake flows show such a similarity in the fully developed state. 

It is convenient to make a distinction between complete and incom- 
plete similarity. In the case of complete similarity , a length 
scale and a velocity scale are sufficient to make the thus-reduced 
velocity distributions coincide completely. In the case of incom - 
plete similarity , more than one velocity scale and/or more than 
one length scale Is required. 

To determine the possibility of the existence of a self- 
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preserving flow, let us set 
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( 23 ) 


Here L q is a characteristic length, and U a characteristic velocity. 

By substitution in the Reynolds equations, it can be shown 
that complete similarity exists only if 

q = 1; p « r « s/2. 

With the relations in (23), it follows from the dynamic equations 
that a plane free jet and a round free .let may have complete simi - 
larity , where n = 1/2 for the plane case, and n = 0 for the round 
case. On the other hand, the plane and the round wake flows can 
have only Incomplete similarity , but in both cases n = 0, and so 
e m is independent of Xj. 

To study memory effects in self-preserving flows, we must 
also determine the possibility of similarity for the relaxation 
equation (16). For this we set 


A ' “ • 
% 


(24) 


Inserting relations (25) and (24) into eqn. (16), we find that 
t = 1 is a necessary condition for similarity. Thus, the memory 
length scale must vary linearly with Xj. For a fully developed 
plane wake flow, s » -1, q 3 1/2 and p « -1/2, and so n = 0 . 


Hence, the memory effect present can no longer be directly 
detected, even if Aj « x, and conditions (22a) and (22b) were not 
satisfied. This follows from-eqn. (10), where for this similar 
flow the effect of the first term can be allowed for with an 
apparently larger e . 

3. Experimental Investigations 

Two different flows with pronounced memory effects are con- 
sidered: the restoration of a strongly disturbed boundary layer 

and the development region of the wake flow behind a cylinder. 

Restoration of a Disturbed Boundary Layer 

The boundary layer in question is that along a straight glass 
plate set up in the plane of symmetry of a wind tunnel. The test 
section of the wind tunnel Is 4.5 m long and has a cross-section 
of 0.8 x 0.7 m 2 . In the experiments the external velocity of the 
boundary layer was kept constant: U Q * 10.5 m/s. A hemisphere 

40 mm in diameter was attached to the glass plate with the midpoint 
of its flat undersurface 3-65 m from the leading edge of the plate. 
The hemisphere was thus entirely submerged in the boundary layer. 
The three-dimensional wake flow of the hemisphere produces a con- 
siderable disturbance of the boundary-layer flow, as shown in 
Fig. 3. 


The velocity measurements were performed with a constant- 
temperature hot-wire anemometer. The heating wire consisted of 
platinum-plated tungsten wire 5 um in diameter. The length of 
the sensing element was 2 mm in most cases, but occasionally 1mm. 
The distance between the wire supports was 10 mm. Measurements 
were made of: the three components of the mean velocity and the 

turbulent velocity, the turbulent shear stresses, spatial velocity 
correlations and one-dimensional spectral distributions of the 
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energy of the axial turbulent component. Wall shear stresses 
were measured with a Preston tube wherever possible. 


Results of these investigations have already been published 
[22]. Here we shall again present some of the most important re- 
sults together with results of further investigations. Equations 
(11) and (16) form the starting point for interpreting the re- 
sults. It was assumed as a second approximation that only the 
change of 3 Uj/3x 2 in the main flow direction [that is A , 3 2 U 1 /3x 1 3x 2 
in eqn. (11)] can be considered. This corresponds to the concept 
that fluid particles which contribute to momentum transfer travel 
only relatively short transverse distances in the x 2 - and x s - 
direction. However, closer investigation has shown that the ne- 
glected terms in eqn. (8) are actually not negligibly small. 

Since the individual contributions have positive and negative 
values, though, the sum total proves to be small. 


An important assumption is now made regarding eqn. (11): 

that the eddy viscosities e and e° are equal and have a constant 

mm 

value in the important x j-direction. 




Fig. 3. Wake flow of a hemi- Fig. H. Distributions of (e^q* 
sphere in a boundary layer. l/tt(c n ) fi and e m acc. to eqn. (il), 

Xj ** 0.50 m, x 2 * 0.02 m. 
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Pig. ** shows e m as a function of x 2 /6 t| for x, ■ 0.5 m and 
x t ■ 0.02 m, first calculated from the Boussirfeaq relation (first 
approximation) 


(«.)■ - - . 


Second, is shown corrected for the influence of intermittent 
m 

turbulence in the outer region of the boundary layer l/ft(e ) 0 , 

m p 

' with ft as the intermittence factor; third, the value is calculated 

from eqn. (11) and is again corrected for intermittence effects. 

Lacking our own measurements of ft, we assumed that the values 

measured by Klebanoff [23] Ju an undisturbed boundary layer were 

approximately valid for the disturbed boundary layer as well. It 

was found that the distribution of e calculated from eqn. (11), 

m 

unlike (e ) fl , is uiaii 1 \r to the distribution in an undisturbed 
m p 

boundary layer-. This L; shown in Pig. 5, which also gives the 
distributionr >.-f the dimensionless quantity e^u^d,, for Xj ■ 

0.25 m, x, ■ ’'.02 in; x 4 * 0.50 m, x s ■ 0.02 m of the disturbed 
boundary layer and for the same * x distances, but x, ■ -0.15 m 
of the undisturbed boundary layer; u* is the shear stress velo- 
city. All values have been corrected with the intermittence fac- 
tor ft, i.e., (-u 2 u x )/ft was used in eqn. (11). The length was 

calculated e« 0 m the equation 




(25) 


The value O.e > f t : e onstant is a mean value of measured 6« 
^lL^lf and a £ rees with the theoretical values found by Saffman 
[2*1]. In addition, A.^ was set equal to 0.*4 in Fig. *4 and 5. 
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/ig. 5. Distributions of the 
dimensionless quantity Cm/u** 

t • calculated from 

(-u 2 u.) measurements of Kle- 
banoff [ 23 ]. 


Fig. 6. Turbulent shear stress, 
measured (+) and calculated by 
Boussinesq relation (o) and by 
eqn. (20) (A). (Boussinesq: 
e m * 0.065 u*6,i ; eqn. (20): 
c m ■ <e m > ■ 1 ^4 . 2 • 10“ ** m 2 /s ) . 


According to our own measurements of the one-dimensional 
spectral energy distribution of u,, i.e. of F.,(kj) (where kj is 
the wavenumber in the Xj-direction) , a value of 0.6 appears to be 
a better approximation than 0.*4. The locally measured value of 
the disturbed boundary layer was taken as u*. Because the quotient 
3 2 U j /3x 1 3x 2 was not constant along a distance equal to A Jt a 
weighted mean was assumed. It was also found that A, was not con- 
stant along a distance equal to A,, as was predicted by eqn. (25). 
Thus, eqn. (16) [or eqn. (20)] is more suitable if one wishes to 
consider a variable Aj. From measurements of U x , uj and fol- 
lows the relation 


Ax K 0,27*^ (A x trod x, in in) . (26) 

Because the result of calculations with various relations for Aj 
has proved to be largely insensitive to changes in this relation, 
the approximate relation shown above is suitable. 

The solution of eqn. (16) with Aj according to eqn. (26) is 
given by eqn. (20) with a ■ 0.27 and n ■ 0.7. Also in eqn. (20), 

x.. « 0 and -u,u,(0) « 0.107 m 2 /s 2 . The calculation was performed 

10 2 1 
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along a line x ■ 0.02 m, x ■ 0 (plane of symmetry). Actually 

the calculation 3hould be carried out along a streamline. But 

because the difference between a streamline and a line x 2 ■ const 

is quite small, the calculation was performed with x 2 ■ const for 

simplicity. According to Fig. 5, e jn /u*6 f| is independent of Xj. 

In eqn. (20) the variable value of u § 6,. should be used for e . 

It was found, however, that for 0.1 m < Xj < 0.6 m, u*6 |f and 

thus e changes only very little, on the order of c « x, 6 * 15 . 
m m * 

Again, for simplicity, a constant mean value ■ 1^.2 « 1 0“ ** m 2 /s 
was used for the calculation by eqn. (20). Fig. 6 shows the re- 
sult c f the calculation of (-u 2 Uj ) together with the measured 
values. The turbulent shear stress was also calculated by the 

Boussinesq formula, with the local value (c )_ * 0.065 . 

m b 

Although the shear stress (-u 2 u ) calculated by eqn. (20) ;.s 
still somewhat smaller than the measured shear stress, eqn. (2 r ‘) 
shows considerable improvement over calculation by the Boussinf-cq 
relation. The shear stress calculated with the local velocity 
gradient 8 Uj/9x 2 is much too small, particularly In the first re- 
gion of the wake flow of the hemisphere. For large Xj, at which 
the disturbed boundary layer is nearly restored, the two calcu- 
lated shear stresses approach the measured value. 

It may be asked whether memory effects in the sense understood 
here might not be present in the undisturbed boundary layer. Re- 
lation (26) does not satisfy condition ( 24 ) with t * 1 for simi- 
larity. Hence, memory effects must be negligibly small. This 
means that conditions (22a) and (22b) must be satisfied. With 
Aj calculated according to eqn. (25), we find for the undisturbed 
boundary layer 
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so that condition (22b) may be satisfied with certaintv. 

Wake Flow of a Cylinder 

This investigation was conducted in a wind tunnel similar to 
that described above. The measurements were made with cylinders 
of various diameters. Here we shall consider only the results ob- 
tained for a cylinder of D = 0.0*1 m. With the velocity of the un- 
disturbed wind tunnel flow T J = 10.5 m/s, the Reynolds number 

Re n = U /D = 26,000. 

Do 

The "classical" investigations of Townsend [8] in a cylinder 
wake (with Re^ = 1360) already showed that self-preservation of 
the mean flow velocity occurs no sooner than Xj/D = 70 to 100, 
because even greater downstream distances are necessary for tur- 
bulence. The investigations performed here pertain to the re- 
gion 2 < Xj/D < 90; the region 20 < Xj/D < 50 was of particular 
interest in terms of the intended memory effects. For Xj/D < 20, 
effects of more or less discrete K£rmcin eddies are noticeable, 
depending on the Reynolds number. 

The hot-wire anemometer described above was used to measure 
the mean velocity components Uj and U 2 , the turbulence intensities 
uj and uj, the shear stress (-u 2 u : ) and spectral distributions of 
Uj. From the measured values of Uj, uj and the spectral distri- 
butions, A 1 (x 1 ,x 2 ) was calculated. 

For x 2 = const, A t = const x x , and thus Aj is a linear func- 
tion of Xj . The eddy viscosity e m was computed for the cross- 
section x 1 = 22.5 D as a function of * 2 /(x 2 ) 99 . Here (x 2 ) 99 is 
the lateral distance for which the velocity difference with re- 
spect to the free stream is one percent of the maximum local velo- 
city difference. The calculation was performed only for the 
inner region, in which the turbulence still shows practically no 
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intermittence. Pic. 7 shows e m /U o D according to Boussinesq, 
according to eqn. ( 11 ) with and according to eqn. ( 1 6 ) . 

The equilibrium value for larger x x , at whiqh the flow is fully 
developed is equal to 0.018. A memory effect is clearly evi- 
dent. The values calculated by eqn. (11) and eqn. (16) are al- 
most equal to this equilibrium value. This value is somewhat 
higher than 0.0164, the value that follows from the measurements 
of Townsend [20] with Re^ = 1360. 
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Fig. 7- Distribution of e / 
U D according to Boussines^ 
(x>, eqn. ( 11 ) (o) and eqn. 
(16) (A). 


Fig. 8 . Distribution of e /U D 
in a cylinder wake along tffe 0 
line x 2 /D = 1.5; calculated ac- 
cording to Boussinesq (o) and 
eqn. ( 11 ) (x). 


The variation of e /U D in an axial direction was calculated 

m o 

according to eqn. (11) along a line x 2 /D = 1 . 5 . This line passes 
through the point in cross-section Xj = 18 D at which the shear 
stress attains a maximum. In eqn. (11) Aj was calculated from 
eqn. ( 25 )» where A^ = 0.5 to 0.6 (x 2 ) 99 according to the spec- 
tral measurements for x 2 = 1.5 D. This leads to A t ^ 0.3 to 
0.36 Xj, and thus the linear relation previously mentioned. Along 
the line x 2 = 1.5 D, the intermittence factor £2=1. Fig. 8 shows 

e /U D calculated according to Boussinesq and according to eqn. 
m o 

(11). After Xj/D = 20 the value calculated from this equation is 
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Tf* (****) 


practically constant and equal to 

the equilibrium value of 0.018. 

The memory effect according to 

eqn. (19) was investigated along 

the same line x 2 * 1.5 D. Here, 

according to Fig. 8, e = const 

m 

along this line. At point x 10 =0, 

it was assumed that -u 2 u x (0) = 0. 

The result of the calculation with 

e = 0.018 U D * 0.0076 m 2 /s and 
m o 

a = 0.3 is shown in Fig. 9, to- 
gether with (-u 2 Uj ) calculated ac- 
cording to Boussinesq with the 

same value of e . This figure 
m 

also shows (-u 2 Uj ) calculated ac- 
cording to eqn. (11) and the 

The linear relation ^ 0.3 Xj can only be an approximation. 
In reality Aj must increase somewhat parabolically with Xj before 
it can vary linearly with x } in the fully-developed equilibrium 
region in which similarity obtains. According to eqn. (25), 

Aj a A lf U Q /uJ, because Uj n. U Q for very large Xj/D. In the simi- 
larity region, A^ f (XjD) 1 / 2 and uj/U o <= (D/Xj ) 1 / 2 ; thus, h 1 « x x . 
From the measurements of Townsend [8] it can be calculated that 
A, ^ 0.26 x r 

With a linear relation for eqn. (16) can satisfy the con- 
ditions for similarity. The following values for the exponents 
in relation ( 23 ) hav •* been obtained for self-preserving wake flow: 
p = -1/2, q = 1/2 and s = -1; thus, n = 0, and e m is independent 
of Xj. As stated earlier, a memory effect is no longer directly 
detectable under these circumstances. To determine the possibility 
of its existence, one must weigh the quantities of the terms in 



Fig. 9. Axial shear stress 
distribution in a cylinder 
wake, measured (•) and cal- 
culated according to Bous- 
sinesq (x), eqn. (19) (o) 
and eqn. (11) (+) . 


measured variation of (-u Uj). 
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eqn. (16) against each other. A rough estimate based on the mea 
surements of Townsend yields 


A x 8 


(“ *4«i) 


ss. 0,2 bi« 0,3 , 


It is apparent that this term is not small. It is possible, then, 
that the memory effect in question is present in the entire wake 
flow, but is not directly detectable due to the similarity of the 
flow. In this case a theoretical solution with the Boussinesq re- 
lation for shear stress and an adjusted value of e without a 

m 

memory effect may be satisfactory. 

Concluding Remarks 

The above theoretical considerations, supported by experimen- 
tal findings, have clearly shown that the additional memory effects 
in turbulence can play an essential role, especially if we are 
dealing with flows which fluctuate relatively rapidly in the flow 
direction, so that the turbulence has no opportunity to adapt com- 
pletely to the local state of the mean flow. Then the Reynolds 
stresses, for example, can be thought of as consisting of two parts 
The first part is practically "local” in character and is deter- 
mined by the local state of the flow, while the second part is an 
additional contribution from the flow's prior history, i.e., an 
additional memory effect due to a relatively rapid change in the 
flow state. These additional memory effects can still be present 
in fully developed turbulent flows which are practically in equi- 
librium. If the flow is self-preserving, however, these memory 
effects are not directly detectable. The above investigations 
suggest that memory effects are always present in self-preserving 
wake flows, because these effects are negligibly small in a fully 
developed turbulent wall flow (boundary layer). 
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When we consider relations which are local in character, 
such as the Doussinesq relation, as a first approximation, the 
above expansions are nothing more t^an a second approximation. 

This point must be reemphasized. 

As stated earlier, memory effects are directly detectable in 
the development region of a turbulent flow, and not only in a 
wake flow, but also in a free jet. An investigation of memory 
effects in the development region of a grid flow is currently 
in progress. It is known that the turbulence of the fully de- 
veloped grid flow is axially symmetrical with u* > u* = u*. If 
these quantities are interpreted as normal stresses, the normal 
stress in the axial direction uj is greater than the normal 
stresses in the transverse direction. Thus, the shear stress 
must be maximum in areas which form an angle tt/ 4 with the axial 
direction. This is incompatible with the absence of a mean velo- 
city gradient, if we assume a simple relation between shear stress 
and deformation rate. The mean flow and turbulence are not yet 
homogenous in the development region, and memory effects can pro- 
duce the aforementioned anomaly, which may exist for an unlimited 
distance downstream provided the memory effects are sustained. 

One difficulty in making a quantitative determination of 
memory effects is our lack of knowledge with respect to relaxa- 
tion time, memory function and Lagrangian autocorrelation. Re- 
lation (25) is nothing more than an approximation. Moreover, the 
determination of the axial derivative of 311^9x2 in eqn. (11) and 
of 9/9Xj (-u 2 Uj ) in eqn. (16) can be performed only with a low de- 
gree of accuracy. The large scatter of the points, as in Fig. 5 
and 7, is due largely to this fact. In this respect a quantita- 
tive determination of the memory effect by means of eqn. (17), 
which requires integration instead of differentiation, offers 
significant advantages. 
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Our considerations here have been limited to turbulent mo- 
mentum transfer. Naturally, such memory effects also play a role 
in the transport of other transferrable quantities such as heat 
and mass. 
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